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Abstract 

The extreme statistics of time signals is studied when the maximum is measured from the initial 
value. In the case of independent, identically distributed (iid) variables, we classify the limiting 
distribution of the maximum according to the properties of the parent distribution from which 
the variables are drawn. Then we turn to correlated periodic Gaussian signals with a 1//" power 
spectrum and study the distribution of the maximum relative height with respect to the initial 
height (MRHi). The exact MRHj distribution is derived for q = (iid variables), a = 2 (random 
walk), a = A (random acceleration), and a = oo (single sinusoidal mode). For other, intermediate 
values of a, the distribution is determined from simulations. We find that the MRHi distribution 
is markedly different from the previously studied distribution of the maximum height relative to 
the average height for all a. The two main distinguishing features of the MRHj distribution are 
the much larger weight for small relative heights and the divergence at zero height for a > 3. We 
also demonstrate that the boundary conditions affect the shape of the distribution by presenting 
exact results for some non-periodic boundary conditions. Finally, we show that, for signals arising 
from time-translationally invariant distributions, the density of near extreme states is the same as 
the MRHi distribution. This is used in developing a scaling theory for the threshold singularities 
of the two distributions. 
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I. INTRODUCTION 

The importance of extreme value statistics 

mm 

has long been recognized in engi- 
neering fields such as hydrology as well as in insurance and finance {g], where gauging the 
effects of catastrophic events is a central concern. Fascination with catastrophic events has 
also brought extreme statistics into the focus of everyday interest, as witnessed by debates 
about climatic events, such as the most violent tornado or the hottest summer of the last 
century [3]. In physics, on the other hand, the use of extreme value statistics has not been 
widespread. The reason for this may be that the rarity of the extreme events naturally puts 
a high price on obtaining information. Nevertheless, the last decade has seen an increas- 



ing interest in extreme value statistics in physical applications, re 



ated, for example, to the 
12, HQ, to fragmen- 



ground state of spin glasses j^, to interface fiuctuations 

tation problems [l5|, to level-density problems of ideal quantum gases [16[, to atmospheric 
physics [17], etc. 

Some of these applications involve extensions of mathematically well known results for 
independent, identically distributed (iid) variables l|, 2, 3, 4] to the physically relevant case of 
strongly correlated variables. This has lead to some exactly solved examples and simulation 
studies of particular systems with well known correlations. A systematic classification of the 
effect of correlations has, however, not yet emerged. 

In the simplest case, extreme value statistics emerges from random numbers drawn from 
a given distribution without any reference to the order of the draws. Often, however, the 
extremum is selected from a well ordered set of random variables. For example, the extremum 
may be the maximum of a time series h{t) in the interval < t < T or, equivalently, the 
maximum height of an interface h{t) in the two-dimensional space {t,h). A relevant point 
to note here is that in correlated systems the boundary conditions at t = and t = T 
may be important. In particular, it has been demonstrated for one-dimensional interfaces 
with periodic or free boundary conditions that, in the presence of strong correlations, the 
boundary conditions do affect the extreme value distribution [l2| . 

The sensitivity to boundary conditions brings up the question whether the extreme statis- 
tics also depends on the zero level from which the maximum is measured. In the iid case, 
the maximum is usually specified with respect to a fixed zero level, related to the scale of 
parent distribution. In other cases, however, one may want to define the zero level through 
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some average measured in the random system. For example, in the case of an interface, 
it is convenient to specify the maximum height of a given reahzation h{t) with respect to 
its average h = h{t)dt, which is also a fluctuating variable 18|. In cases where the 

fluctuations of h diverge in the limit T — oo, it is not surprising that the extreme height 
distribution is sensitive to the choice of the zero level. 

In physical applications, choosing the origin at the first value of the measurement, i.e., 
measuring the maximum height with respect to the initial height, is one of several natural 
possibilities. In finance, as well, one may be interested in the probable extremes of stock 
prices with respect to the starting price. Measuring the height of a signal from the initial 
height instead of the average height, for example, may seem like a trivial shift of the origin. 
However, if the initial height is a random variable, with a distribution determined either by 
the experimental setup or by the inherent statistical properties of the infinite signal, then 
the extreme statistics may depend on the probability distribution of the initial value. 

In this paper we study the distribution of the quantity m = maxt[/i(t) — h{0)], i.e., of the 
maximum height relative to the initial height (MRHi). After a comprehensive review of the 
iid case, we investigate m for correlated Gaussian signals h{t) with a 1//" power spectrum 
and with periodicity h{t) = h(t + T). It is interesting to compare the MRHj distribution 
with the distribution of the quantity hm = maxf[/i(t) — /i], i.e., of the maximum height 
relative to the average height (MRHa), recently analyzed in Refs. Q,Q- One of our main 
findings is that the MRHj and MRHa distributions are different for all a. 

The paper is organized as follows: The case of independent, identically distributed (iid) 
variables is considered in Sec. II. We show that the MRHj distribution is the same as the 
recently studied density of near extreme events 20], and we find that the MRHi and MRHa 
distributions may or may not differ, depending on the tail of the parent distribution. The 
model of periodic, correlated, Gaussian signals h{t) with a 1//°^ noise spectrum is introduced 
in Sec. HI, and our procedure for determining the MRHi distribution from simulations is 
described. In Sec. IV the exact MRHi distribution is derived in the special cases a = (iid 
Gaussian variables), a = 2 (random walk), a = 4 (random acceleration), and a = oo (single 
sinusoidal mode), including the dependence of the distribution on the boundary conditions 
for a = 2 and a = 4. Simulation results for the MRHi distribution with periodic boundary 
conditions for other, intermediate values of a are presented in Sec. V, and the evolution 
of the distribution with changing a is discussed. In Sec. VI we show that the equivalence 
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of the density of near extreme events and the MRHi distribution, demonstrated for iid 
variables in Sec. II, continues to hold for correlated signals, which are time-translationally 
invariant. This is then used to determine the scaling of the singularities of the two equivalent 
distributions. Finally, Sec. VI contains concluding remarks. 



II. MRHi DISTRIBUTION FOR INDEPENDENT, IDENTICALLY DIS- 
TRIBUTED VARIABLES 

Consider N random variables hi, h2, h^, selected independently according to the 
probability density p{h), called the parent density. The probability that the variable hi is 
less than x is fi{x) = J^^p{h) dh, so the probability M^g,x{x, N) that all N variables take 
values less than x is 

M^a.(x,iV)=/i^(x). (1) 

Since Mmax(a;, A^) also represents the probability that /imax = maxj hi is smaller than x, the 
probability density of the maximum is 

Pmax(a;, N) = -^ M^ax(x, N) . (2) 

The asymptotic form of the distribution function ([2]) for large is discussed in stan- 
dard textbooks 3|, |j] on extreme value statistics. For a wide class of parent distributions, 
Mmax(2;, N) becomes independent of in the large N limit, on making the linear change of 
variable x = a^z + b^^ with suitable parameters a^, b^: 

M^^^{a^z + b,,N)^ M^^^iz) , P;,,(^) = ■ (3) 

Depending on the tail of the parent distribution, the limit function M^^^{z) belongs to one 



of three classes, associated with the names of Fisher-Tippett-Gumbel, Frechet, and WeibuU. 



We now turn to the main subject of this paper, the statistics of the quantity /imax — ^i, i-e., 
of the maximum height relative to the initial height. For identically distributed variables 
it does not matter which hi is singled out as a reference, but to be specific we use hi. 
Let us calculate the probability J^{m,N) that all of the relative variables hi — hi, /12 — 
hi, . . . , hN — hi are less than m. Since the first of these relative variables is identically zero, 
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jF(m, A^) vanishes for negative m and for positive m is the same as the probabihty that 
h2 — hi, . . . , Hn — hi are all less than m. Thus, 

dhip{hi) / dh2p{h2) . . . / dhj^p{h^) , 

oo J ~oo J ~oo 

/oo 
dhi p(/ii)M^ax(/ii + m, AT - 1) , (4) 
-oo 

where Q{x) is the standard Heaviside function, and in going from the first line to the second 
we have used Eq. ([T]). Differentiating J^{m,N) in Eq. (jlj) with respect to m and making 
use of Eq. ([2]) we obtain 

1 

P{m, N) = - 5{m) + e(m) / dhpih) P^,^{h + m, iV - 1) , (5) 

' J — oo 

for the probability density of the maximum of all relative variables. 

To analyze the asymptotic form of P(m, A^) for large A^, we neglect the first term on the 
right-hand side of Eq. ([5]) and substitute the hmiting distribution 

P^ax(x, AT - 1) ^ a-'P:,,^{a-\x - h,)) , (6) 

introduced in Eq. ([3]), in the second term. The convolution of the two normalized distribu- 
tion functions P* and p in Eq. ([5]) depends on their relative scales, and we distinguish the 
three cases: (i) vanishes, (ii) converges to a finite a value, and (iii) diverges. In the 
first case Pmax(a;, A^ — 1) ~ 5{x — b^), whence P(m, A^) ^ p{m — hj^). In case (ii), Pmax and 
p vary on the same scale, and the MRHi distribution is a convolution of two non-degenerate 
functions. As we shall see later, in this case P^axl-^) Fisher- Tippett-Gumbel form, 

P^^yr{z) = exp {—z — e^^). Finally, in case (iii), aj^p{aj^z) ~ S{z), and the integral is readily 
evaluated on changing the integration variable to z. The results can be summarized as 



P(m, A^) 



R{m-h^), a^-^a (7) 



oo 



where 



R{x) 



/oo 
dzp{az — x) exp (— z — e^^) . (8) 
-oo 



We now review the connection between the sea 
important role here, and the parent distribution 



e parameters a^, b^, which play such an 



It is useful to work with the 



functions g{x) and f{z), defined by 

/^(x) = exp [-e-^(-)] , M;;,,(z) = exp [-e-^(^)] . (9) 

Here is the integrated parent distribution introduced above Eq. ([1]), and Mj^^^(z) is the 
hmit function defined in Eq. ([3]). Note that g diverges as /i approaches 1, and its asymptotic 
form is related to the tail of the parent distribution. According to Eqs. ([3]) and Qj, 

Cax(^) = ^ exp [-e-^(^)] . (10) 

In terms of g{x) and f{z), the A^-independence (|3]) of the extreme distribution for large 
N takes the form 

g{a^z + bJ-\nN^f{z) . (11) 
We define the scale factors and by 

gib,)=\nN, g\b,) = a'^\ (12) 

Here the conditions /(O) = and /'(O) = 1 have been imposed for convenience, but may be 
relaxed by making a linear transformation z = az' + /?, with parameters a and /3 which are 
independent of A^. Such transformations do not change the asymptotic behavior of and 

^- 

We now relate the three cases in Eq. ([7]) to the tail of the parent distribution. According 
to dn]), 

which implies that a^r (i) goes to zero, (ii) converges to a constant, or (iii) diverges if the 
asymptotic dependence of ^at on In A^ is slower than linear, linear, or, faster than linear. From 
Eq. ( fT2l) we see that in these three cases g[x) diverges, for x — > oo, (i) faster than linearly, 
(ii) linearly, or, (iii) slower than linearly, and since /i(x) ^ 1 — e^^'^^^ for large g, 1 — ^{x) 
vanishes (i) faster than exponentially, (ii) exponentially, or (iii) slower than exponentially, 
respectively. Here "exponentially" is understood in the broad sense as corresponding to a 
linear leading divergence of g{x). Case (ii) therefore includes parent distributions for which 
1 — fi{x) decays as x^e~'^^ or e*^' e~^^, e < 1, since, for all these cases, g{x) ~ cx to leading 
order as x —>■ oo. 

Finally, recalling that p{x) = dfi{x)/dx, we conclude that the MRHj distribution is given 
by the second line on the right-hand side of Eq. ([7]) whenever the parent density distribution 



p{x) decays exponentially for large x, in the broad sense just defined. In this case -Pmax is of 
Fisher-Tippett-Gumbel form as used in Eq. [71 If the decay is more rapid than exponential 
or less rapid than exponential, then lines one and three in Eq. ([7]) apply. 

It is enlightening to supplement this general discussion with explicit results for some 
characteristic parent distributions. We begin with the case of a parent distribution p{x) 
which for large x decays according to the generalized exponential form x^exp [— (x/^)*^], 
where 6 > but A can have either sign. According to Eq. (Q, g{x) has the asymptotic 
form g{x) ~ {x/^Y for large x to leading order, independent of A, and from Eqs. (fTTj) and 

f{z) = z , — oo < 2; < oo , (14) 
a, ^ ib-\\^N)^l'-^ , ^ e(lniV)i/^ . (15) 

From Eq. f|T5|) . we see that 0, a, and oo for 5 > 1, 5 = 0, and 5 < 1, respectively. 

Thus, the corresponding MRHi distributions are given by the first, second, and third lines, 
respectively, on the right side Eq. ([7]), in agreement with the general conclusions of the 
preceding paragraph. Substituting the scaling function (fl^ in Eq. (ITOl) yields the Fisher- 
Tippett-Gumbel form of the extreme distribution -Prnax(^)' already shown explicitly just 
above Eq. ([7]) and in Eq. ([8]). 

Next we consider a parent p(x) which vanishes for x greater than a finite value X\ and 
varies as {x\—x)°'^^, with a; > 0, as x approaches X\ from below, so that 1 — /i(x) ~ (xi— x)", 
and g[x) ~ — aln(xi — x) for small x\ — x. From Eqs. (|TT1) and (|T2l) . 

f{z) = -a\\\{\-a~^z) , < 2; < a , (16) 
^ a-^N-^/" , b^^xi- N-^'"" . (17) 

Since a^v — > in the limit — 00, the MRHj distribution is given by the top line on the 
right side Eq. ([7]). This is consistent with the general analysis given above, since the parent 
distribution has a faster than exponential decay for x ^ 00, having already attained at 
the finite value xi. Substituting the scahng function ( |T6l) in Eq. ( ITOj) yields the WeibuU 
form of the extreme distribution -Pmax(^)- 

Finally we consider a parent distribution p(x) which for large x decays as x~^~^, with 
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P > 0, so that 1 — fi{x) ~ X ^. In this case g{x) ~ p\nx, and from Eqs. (fTTi) and (fT2i) . 

/(z) =/31n(l + /3-^z) , -/3<z<oo, (18) 
^ p-^N^/'^ , ^ N^/f^ , (19) 

Since a^v — > 0, the MRHi distribution is given by the third hne on the right side Eq. ([7]). 
This is also consistent with the general analysis given above, since the parent distribution 
has a slower than exponential decay. Substituting the scaling function (fTSll in Eq. (fTOj) 
yields the Frechet extreme distribution. 

The expressions for f{z) in Eqs. (fT4l) . (fT6l) . and (fTSl) . correspond to the special cases 



77 = 0, ?7 < and rj > 0, respectively, of the generalized extreme value distribution 2l| with 
f{z) = ?7ln(l +r]~^z). Here we have treated the three cases separately in order to highlight 
the different asymptotes of qn and Bn- 

Recently, the density of near extreme states, corresponding to the distribution of the 
relative variables /imax — hi and defined by 

pir,N) = - J2 Hr-iKn..-h)] , (20) 

was studied by Sabhapandit and Majumdar for iid variables. Since the MRHi distribu- 
tion P(m, A^) = {6 {[m — (/imax — ^1)]) that we consider does not depend on the particular 
variable hi chosen as a reference, it may be rewritten as 

1 ^ 

P(m, N) = -J26[m- {hm.. - hi)] = N-'5im) + p{m, N) . (21) 

i=l 

Thus, the density of near extreme states p{m, N) and the MRHi distribution P{m, N) 
only differ by a delta function. The results of Sabhapandit and Majumdar for the limiting 
behavior of p(m, A^) for large are essentially the same as ours for P{m,N), apart from 
our more general evaluation of the threshold case of exponentially decaying parent. 

In Section VI we point out that the equivalence between the density of near extreme states 
and the MRHi distribution is not limited to iid variables, but also holds for the periodic, 
correlated signals considered in Sec. HI- VI. 
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III. CORRELATED GAUSSIAN SIGNALS 



We now turn to Gaussian signals of periodicity h{t) = h{t + T) with configurational 



weight 



19 



22 



V[h{t)] oc e-^l'^Wl , (22) 
where the effective action, in Fourier space, is 

N/2 

5[cfc;a] = (27r)-Ti-"^n"|c„p. (23) 

n=l 

Here the Ck are coefficients in the finite Fourier series 

N/2 

h{t)= c„e2™*^^' < = c_„, (24) 

n=-N/2+l 

where is a positive, even integer. Since the maximum frequency appearing in the sum is 
of order N/T, the series does not resolve fine structure on a time scale less than r = T /N . 

We will be mainly interested in the continuum limit — > oo, r — with T = Nt fixed. 
Expressed in terms of h{t) instead of its Fourier transform, the action in Eqs. ( I22l) and ( 123|1 
takes the form 



S[h{t)] = \j^ dt 



in this limit, which implies the stochastic equation of motion 



(25) 



mm) = s{t-t'), (26) 



where ^(t) is Gaussian white noise with zero mean. 

The requirement c„* = c_„ in Eq. (!24l) guarantees that h{t) is real. The Fourier co- 
efficients Co (present for even only) and c^^^ real, but the other c„ are complex. 
Configurational averages involve integration with the statistical weight (122|) . (123!) over the 
phase space 

/oo N/2— 1 ^ ^ 

n / dRe[ck] / dlm[ck] . (27) 
-oo k=l J —oo 



From Eqs. ( I22l) and (l23l) one sees that the amplitudes of the Fourier modes are inde- 
pendent, Gaussian distributed variables, but only for a = are they identically distributed. 
This is also apparent from the mean square amplitude (|c„p) oc which is consistent 
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with a 1//" power spectrum and independent of n only for a = 0. Tuning a allows us to 
treat a broad range of time signals and recover some important special cases. The values 
a = 0, 1, 2, 4 correspond, respectively, to white- noise (iid Gaussian variables), 1// noise |23|, 
the random walk (diffusion), and the random acceleration process 2J]. For a = 0,2,4 this 
correspondence is immediately apparent from the stochastic equation of motion (1261) . 

Although the Fourier components c„ are uncorrelated, the corresponding time signal h{t) 
is correlated at different times t and t' for a > 0, and the correlation increases with increasing 
a. For example, for a = the h{t) at different times are iid random variables, whereas for 
a — i> oo, h{t) becomes a single- mode sinusoidal curve, as discussed below. For < a < 1, 
the correlation function {h{t')h{t' + t)) is bounded, while for a > 1 it diverges in the limit 
T — >■ oo with t/T finite. For a more detailed discussion of correlations in 1/f" signals, see 

Q. 

In the next two Sections we study the distribution function of m = maxt[h{t) — h{0)], 
i.e., of the maximum height with respect to the initial height (MRHi), for the Gaussian 
model defined by Eqs. ( |22l) and ( |23l) . First we derive the MRHi distribution exactly in the 
special cases a = 0, 2, 4, and oo and then, for other, intermediate values of a determine the 
distribution with numerical simulations. 

In our simulations the distribution P{m, N) of m was computed from about 10^ to 10^ 
independent signals h{t) for each of the various values of and a that were considered. 
Each signal was generated by selecting the Fourier coefficients c_^^2+i • • • 5 randomly 
from the Gaussian distribution ( l22i) . ( I23l) and then summing the series ( l24l) to obtain h(t) at 
t = 0, T/N, 2T/N, , . . . , (A^ — 1)T/N. For each signal the maximum m of these A^ heights 
relative to the initial height h{0) was determined, and then the values of m for the 10® to 10^ 
independent signals were binned to obtain the distribution P{m,N). This procedure was 
carried out for increasingly large values of A^, to obtain the best estimate of the distribution 
in the continuous time limit N ^ oo, t —>■ with T = Nt fixed. 

To extract scaling functions ^{x), free of fitting parameters, from P{m,N) in the limit 
A^ — > oo, we follow the same procedure as Gyorgyi et al. [13]. In cases where (m)^ = 

mP{m,N)dm and = {{m — {m-) ^)'^)^J'^ have the same large A^ behavior, i.e. where 
their ratio approaches a constant for large A^, we scale by the average, introducing the 
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variable 



and defining the rescaled distribution by 

$(a;) = lim {m) ^P{{m) ^^x, N) . (29) 

Since m = ma.xt[h{t) — h{0)] is non-negative, the distribution function $(a;) is only defined 
for X > 0. According to Eqs. ( l28l) and ( l29l) . is normahzed so that dx^{x) = 1 and 
has the mean value (x) = J^^ dx x^{x) = 1. 

If, on the other hand, (m)^ and 0"^^ have different large behavior, we scale by the 
standard deviation {a scaling), introducing the variable 

y = ^I^^lM^ (30) 

and defining the corresponding distribution function by 

^{y) = lim a^P{{m)^ + a^y,N) . (31) 

According to Eqs. fl30l) and fl3T]) . <!>(?/) is normalized so that J^^dy ^{y) = 1 and has the 
moments {y) = J^^dyy^y) = 0, {{y - {y)fy^'^ = 1. 

IV. EXACT RESULTS 
A. Special case a = 

As mentioned below Eq. fl27j) . the correlated Gaussian signals h{t) defined in the preced- 
ing Section reduce to iid variables in the limit a 0. Since the Gaussian parent distribution 
has a faster than than exponential decay, the MRHi distribution, given by the upper entry 
in Eq. ([7]), is also Gaussian. Since {m)^ ~ (InA^)^/^ and ~ (lnA^)~^/^ for large A^, we 
specify the distribution in terms of the variable y and the function $(?/) of a scaling, defined 
in Eqs. (I30l) and (|3T1): 

l>i(y) = (27r)-i/2 exp(-2/V2) , -oo < < oo . (32) 



In contrast to Eq. ( I32l) . the limiting MRHa distribution, considered in Ref. 1^, has the 
Fisher-Tippett-Gumbel form given by Eqs. (fTOl) and (fT4|) and shown explicitly just above 
Eq. ([7]). The two distributions are compared in FigiH 
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FIG. 1: Results for the MRHj distribution from simulations, plotted using cj-scaling (Eqs. (|30p and 
(j3ip ). for Q in the range < a < 1. The simulations were performed with = 4096 terms in 
the Fourier series (j24p . The exact analytical result (j32p for N ^ oo (normalized Gaussian) is 
also shown, and the simulation results are within the linewidth of the Gaussian. The curve FTG 
corresponds to the Fisher-Tippett-Gumbel distribution, which is the limiting MRHa distribution 
in the interval < a < 1. 

The simulation results shown in Fig. [1] will be discussed in the next Section. It turns out 
that the MRHj distribution has the Gaussian form (1521) not just for a = 0, but throughout 
the interval < a < 1. 



a =0.0 




B. Special case a = 2 

As mentioned after Eq. fl27|) . the correlated Gaussian signals defined in the preceding 
Section correspond, for a = 2, N oo, to random walks h{t) governed by the stochastic 
equation of motion fl2B]) . The statistical weight or propagator Z{h, hQ,t) for a random walk 
from initial position to /i in a time t satisfies the diffusion equation {d/dt — D d^/dh"^) Z = 
0, with initial condition Z{h, ho, 0) = 6{h — hQ). We set -D = 1/2 from now on, corresponding 
to Gaussian signals normalized as in the preceding Section. The particular value of D is not 
important, as it drops out of the MRHi distribution on scaling by the average, as in Eqs. 
(125]) and (I2U]). 

We will need the well known solutions of the diffusion equation 

Mh, ho, t) = j^^yj^ exp [-{h - hof/2t] (33) 
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for random walks in the unbounded space —oo < h < oo, and 



Zi{h, ho, t) = Zo{h, ho, t) - Zo{h, -ho, t) (34) 



for random walks in the half space h > 0, with an absorbing boundary 25| at h = 0. 

Let us consider the family of random walks h{t) in the unbounded space —oo < h < oo 
which satisfy the periodic boundary condition h{0) = h{T). With no loss of generality in 
the result for the MRHj distribution, one may choose /i(0) = 0. The fraction jFpei.(m, T) 
of the walks with endpoints /i(0) = h{T) = which never exceed height m in the interval 
< t < T may be expressed as 

^ Zi{m,m,T) ^ ^ _ 2m2/T (OK) 

Zo(0,0,T) ■ ^ ^ 

Here Zo and Zi are the whole and half-space propagators of Eqs. fl5^ and ([31]), and 
Zh<m{hi, ho, T) is the propagator for random walks from ho to hi in time T which lie entirely 
in the subspace h < m. In going from the first line of Eq. (!35|) to the second, we have used 
the relation Zh^m{hi, ho, T) = Zi{m — hi,m — ho, T), which follows from invariance of the 
statistical weight under the coordinate transformation h m — h. 

Differentiating jFper(m, T) in Eq. ( l35l) with respect to m yields the MRHj distribution 

PMm,T) = ^e-^-'/^ (36) 

for random walks with periodic boundary conditions. The mean value of the MRHj dis- 
tribution (l36l) is (m) = {nT/SY^"^. In terms of the variables x = m/{m) and $(x) = 
{m)Pper{^,T) of Eqs. ( 128|) and ( 129|) . the distribution takes the form 

$i(a;) = ^xe-""'/^ (37) 

The mean value (m) of the MRHj distribution given in the preceding paragraph is the 
same as the mean value (/im) of the MRHa distribution for a = 2 obtained by Majumdar 
and Comtet ^]. The equality 

(m) = {h^) (38) 

is a general consequence of the definitions m = max^ \h{t) — /i(0)] , hm = max^ [h{t) — h] , where 
h is the time average of h{t), and time-translational invariance, on averaging over all paths, 
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in the form {h{0)) = {h). Note that the equahty fl38|) holds for all a, not just a = 2, and for 
non-periodic as well as periodic boundary conditions, as long as the boundary conditions 
are consistent with the time translational invariance. 
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FIG. 2: The exact MRHj distribution (j37p for a = 2 and periodic boundary conditions, and the 
corresponding MRHa or Airy distribution [2]. Both distributions are scaled with the average, as 
in Eqs. ([28]) and ([291). 



As in the case a = 0, the MRHj and MRHa distributions for a = 2 are not the same. 
Majumdar and Comtet 12| have shown that for random walks with periodic boundary 



conditions, the MRHa distribution is the so-called Airy distribution. The two distributions 
are compared in Fig. [2l We see that, for both small and large x, $i(x) has the greater weight. 
These features are found for all a > 1 and can be understood heuristically as follows: 

In the case of ^a{x), small x means small hm = /^max — h, or /imax ~ h. There are very 
few such configurations, only those for which h{t) is nearly constant. In the case of 
on the other hand, small x means small m = /imax " ^(0), or /imax ~ h{0). This condition 
is far less restrictive, since it is satisfied by any configuration h{t) which is below h{0) most 
of the time. Thus, small m is much more probable than small hm, i.e., ^ ^a{x) for 

small X. 

Turning to the large x behavior, we note that for a configuration h{t) which makes a 
very large positive excursion with respect to the initial height h{0), the average height h 
also tends to be much larger than h{0). Thus, for such a configuration m ^ h^- Roughly 
speaking, this means that a large value of hm has the same probability as a much larger value 
of m. This, together with a probability distribution that decreases rapidly with increasing 
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m implies $i(x) ^ $a(3;) for large x, as seen in Fig. [2j 

The distribution of the maximum height not only depends on the reference height from 
which it is measured. It also depends on the boundary conditions imposed on h{t) at t = 
and t = T. Before leaving the special case a = 2, we derive the MRHi distribution for 
random walks for two non-periodic boundary conditions of general interest. 

Consider the family of random walks h{t) on the infinite interval — oo < /i < oo with 
fixed endpoints h{0) = and h{T) = hi. For this fixed boundary condition Eqs. fl35|) and 
(!36|) are replaced by 

{hi,0,T) 



J^{i^{m,hi;T) 



and 



Zo{hi,0,T) 
Zi(m — hi,m, T) 
Zoihi,0,T) 

2(2m - hi] 



1 — exp 



hj - (2m - hif 
2T 



exp 



hj - (2m -/ii)^ 
2T 



(39) 



(40) 



which reduces to Eq. (l36l) for hi = 0. 

Finally, we consider the family of random walks h{t), < t < T, with initial condition 
h = a.t t = but with no restrictions on h{T). For this boundary condition Eqs. (|35ll and 
( l36l) are replaced by 

JTo^dhi Zh<rnihi,0,T) 



JZ,dhiZo{hi,0,T) 
dh Zi{h, m, T) = erf 



/m 
dhi Zi{m — hi, m, T) 
-oo 



m 



(41) 
(42) 



where erf denotes the error function 26||. The quantity jFfi.ee(m, T) is the probability that a 
random walk which begins at the origin has not yet reached point m after a time T. The 
integral on the right-hand side of Eq. (H2!) is the "persistence" probability that a random 
walk which begins at point m has not yet reached the origin after a time T. Equation fH2|) 
states the obvious fact that these two probabilities are equal and reproduces the well known 
T~^/^ decay of the persistence probability for long times. 

Differentiating jFfi.ee(m, T) in (H2l) with respect to m yields the MRHj distribution 



Pfree{m,T) 



2 \ 



1/2 



: 



-m?/2T 



ttT J 



(43) 
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After scaling by the average, as in Eqs. ([2SD and dSHD, with (m) = (2T/7r)V2^ the distribution 
function (l43i) takes the form 

= -e-"'/^ (44) 

TT 

which is clearly different from the result (137|) for periodic boundary conditions. 



C. Special case a = 4 

As pointed out below Eq. fl27j) . the correlated Gaussian signal hit) of Sec. Ill may 
be interpreted, for a = 4, — cxo, as the position of a particle which is randomly ac- 
celerated according to the stochastic equation of motion fl26l) . Following the same ap- 
proach as in the preceding Subsection, we work with the statistical weight or propaga- 
tor Z{h,v] ho,vo]t) for a randomly-accelerated particle with position and velocity hQ,VQ at 
t = and values /i, f at a later time t. This quantity satisfies the Fokker-Planck equation 
[d/dt + vd/dh — Vd'^/dv^) Z = 0, which is basically a diffusion equation for the velocity, 
with initial condition Z{h,v; h(),vo; 0) = 6{h — ho)6{v — vq). We set "D = 1 from now on. 



as in Ref. [2j]. The particular value of V is not important, as it drops out of the MRHi 
distribution on scaling by the average, as in Eqs. ( 128|) and ( l29l) . 

In calculating the MRHi distribution, we will need the solutions Zo and Zi to the Fokker- 
Planck equation in the unbounde d sp ace — oo < /i < oo and in the half space h > 0, with 



an absorbing boundary condition [25|, respectively. Expressions for Zo{h, t>; ho, vq; t) and for 
the Laplace transform Zi{h,v; ho,vo; s) = dte~'^^Zi{h,v;ho,vo;t) are given in Ref. [24| 
and in Appendix A of this paper. 

Let us consider the family of trajectories h{t) in the unbounded space — oo < h < oo with 
periodicity h(t) = h(t + T). As in the case a = 2, we may choose h{0) = 0, with no loss 
of generality in the result for the MRHi distribution. The velocities at t = and T are the 
same, as follows from the periodicity, but otherwise unrestricted. All values of the initial 
velocity are assumed to be equally probable. For this boundary condition, Eq. fl35l) in our 
treatment of the random walk is replaced by 



J^^dv Zi{m,v;m,v;T) 
/^^^^o(0,w;0,f;T) ' 



(45) 
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Here Zh<m{hi,vi; ho,vo',T) is the statistical weight for a randomly accelerated particle that 
propagates from point hQ,vo in phase space to hi,vi in a time T without leaving the half 
space h < m. In going from the first line of Eq. (1451) to the second, we have used the relation 
Zh<m{hi,Vi] ho, Vq] T) = Zi{m — hi, — f i; m — Hq, —Vq; T), which follows from invariance of 
the statistical weight under the coordinate transformation h ^ m — h. The probability 
distribution of m is obtained by differentiating jFpei.(m, T) with respect to m. 



Using the results for Zq and Zi in Ref. 



24| . we have derived the MRHj distribution for 



the periodic boundary condition of the preceding paragraph. The calculation is outlined in 
Appendix A. For the mean value, one obtains (m) = ^ (vrT'^)^''^. In terms of the scaling 
variable x and scaling function <l>(x) of Eqs. (|28l) and ( l29l) . the distribution is given by 

$i(x) = 2-^/^3^^-^^ X-'/' exp (-^ x') ^ (- J4 ' ^^^^ 



where U{a,b,z) is Kummer's function 26|. The function $i(x) shown in Fig. [3] has the 



asymptotic forms 26 1 



2-4/33l/6^-2/3r(|)^-l/3 ^ X^O, 



$i(x) ^ r ^ (47) 
i exp (- ) , X — > oo , 



and the moments 



(-^) = lf^r(|)-% (48) 

for arbitrary u > —2/3. 

The difference between $i(x) and ^a{x) in Fig. [3] is even more dramatic than for a = 2. 
For a = 4, $i(x) diverges at x = and decreases monotonically with increasing x. The 
increasing weight, with increasing a, of the MRHj distribution $i(x) will be discussed below, 
in connection with simulation results for a broad range of a. 

Next, we consider the MRHj distribution for a randomly accelerated particle with position 
h = and velocity Vq, but with no restrictions on the position and velocity at t = T. For 
this boundary condition Eqs. (14T!) and (142|) are replaced by 

^ . ^^ J'"^dhiJ^^dviZh^rn{huVi;0,Vo;T) 

J^tee{m,vo;T) = — — -— (49) 

J_oo dhi dvi Zo{hi,vi] 0, Vq] T) 

/m POO 
dhi I dvi Zi{m — hi, ~vi;m, —Vo;T) (50) 
-oo J — oo 

/»oo /»oo 

= / (i/i / dv Zi{h,v;m, —V(j;T) = q{m, —Vq;T) . (51) 

Jo J-oo 
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FIG. 3: The exact MRHt distribution (1461) for a = 4 and periodic boundary conditions, and the 

□ 

corresponding MRHa distribution, determined from simulations [19|]. Both distributions are scaled 
with the average, as in Eqs. (j28p and (|29p . The asymptotic forms of the MRHj distribution for 
small and large x are given in Eq. (j47|) . The MRHa distribution is exponentially small for small 
X and for large x is essentially Gaussian with a smaller variance [2]. 

Here jFfi.ee(m, fo; T) is the probability that a randomly accelerated particle which begins at 
the origin with velocity vq has not yet reached point m in a time T. As in the case a = 2, 
J-'free represents a "persistence" probability. The quantity q{m, —vq; T) on the right-hand 
side of Eq. f lM]) is the probability that a randomly accelerated particle with initial position 
m and initial velocity —Vq has not yet reached the origin after a time T. Equation flSip 
states the obvious fact that these two probabilities are equal and reproduces (see Eq. (!54|) 
below) the well-known T~^/^ decay 2J] of the persistence probability for long times. 

Combining P^-ceeijn^VQ^T) = di^Qe^{m,VQ]T) / dm, Eg. (ISTl) . and the expression for the 
Laplace transform of g(m, — t>o; T) in Eq. (19) of Ref. |2J], we obtain 



Pfree(m,t;o; s) = / dTe ^'^ P{m,vo;T) 
'o 

f oo 



poo 

/ dF F-2/3e-^-Ai {-F'/^vo + F-'/'s) 
Jo 



X 



47ri/2 



r (-1, |F-.3/.) 



(52) 



(53) 



for the Laplace transform of the MRHj distribution. Here Ai{z) and T{a, z) are the standard 
Airy and incomplete Gamma function. In principle, Pfi.ee(^, ^"0, T) can be determined from 
Eq. (|53|) by integrating over F and inverting the Laplace transform numerically. For = 
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all of the moments (rrf) of the distribution can be calculated analytically 28| from Eq. (l53l) . 
Equation (21) of Ref. 2J] leads to the exact asymptotic form 



free 



;m,fo;T) 



2-l34/3^-3/2p^l) 
' 2\ 1/12 

^) ^ 



d_ 

dm 



1 2 

6'3' 



9m 



(54) 



for T ^ m?^^ and T ^ |fo|^''^. Here t/(a, 6, z) is Kummer's function 26|. 

Finally, we note that the MRHi distribution for a = 4 and all its moments can be 
calculated analytically [28^ for two additional boundary conditions: In the first case h(T) = 
h{0) and v{T) = v{0) = 0. In the second case h{T) = h{0), while v{T) and v{0) are 
independent and unrestricted, and integrated over all values between —00 and 00. 



D. Special case a = 00 

According to Eqs. fl221) - fl2^ . all but the modes with Fourier coefficients c_i and Ci are 
suppressed in the limit a — cxo, so that 

h{t) - h{0) = ci (e^-**/^ - 1) + cl (e-2-^*/^ - 1) (55) 
= 2|ci I [cos (27rt/T + (/?) — cos*/?] , (56) 

where Ci = Icile*"^. Thus, m = ma.Xt[h{t) — h{0)] = 2|ci|(l — cosip) = 4|ci| sin^((y9/2), and, in 
accordance with Eqs. (!22|) and (123|) . the distribution of m is given by 

\ /'27r poo ^ . 

P{m,N) = - dip d|ci| |ci|e-^l'=^l'5 (^m-4|ci|sin2|j , (57) 

^^r,,ln.,re.J--^). (58) 

oTT Jq \ 16 sm V J 

Here we have absorbed several constants in A, rewritten the integral 
(iRe [ci] (ilm [ci] over phase space (see Eq. (1271) ) in terms of polar coordi- 
nates (|ci|, ip), performed the integration over |ci|, and made the substitution = ^ip. By 
carrying out the m integration before the integration, it is easy to check that P(m, N) 
in Eq. ( l58l) satisfies the normalization condition dm P{m, N) = 1 and has the average 
value (m) = (tt/A)^/^. 

The integral over 1!} in Eq. (l58l) can be evaluated with the help of Ref. or Mathematica. 
In terms of the variable x = m/ (m) = (A/7r)^/^m and scaling function $(a;) of Eqs. (|28|) and 
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(129|) . the MRHi distribution for periodic boundary conditions and a = oo takes the form 



2\ -1/4 / 2\ /II 



where again U{a,b,z) is Kummer's function [26]. The function has the asymptotic 



forms 



26| 



and the moments 



$I(X) ^ <( _ ' (gQ) 

2"^/^ exp (— Y^x^) , X ^ oo , 

(^")- 33^,1)-" (61) 

for arbitrary z/ > — |. 

From Eq. (1601) we see that MRHj distribution function (I59p diverg GS cLS X ' for X — s> 0, 
even more strongly than in the case a = 4 (see Eq. (H71) ). where there is an x~^^^ divergence. 
The distribution function fl59l) is plotted in Fig. HI The MRHa distribution for a = oo, 
calculated by Gyorgyi et al. 19(] , has precisely the same form fl37|) as the MRHi distribution 
for a = 2. Thus, the curves a = 2 and oo in Fig. H] correspond to the MRHa and MRHj 
distributions, respectively, for a = oo. 

V. RESULTS OF SIMULATIONS 

In addition to the analytical results for the special values of a = 0, 2, 4 and oo, we have 
determined the MRHj distribution for intermediate values of a by numerical simulations, 
following the procedure described in Sec. Illli This allows us to address some interesting 
questions concerning the a dependence. In the simulations the number of terms in the 
Fourier series (HM was chosen large enough so that finite-size effects were negligible, and the 
exact results of the preceding Section, corresponding to = oo, could be reproduced with 
"width-of-the-line" accuracy everywhere except in the immediate neighborhood of x = 0, 
where $i(x) diverges for a greater than a critical value close to 3. 

First we turn to the weak-correlation regime < a < 1, where the correlation function 
|7i^0)/i(t)) remains finite in the limit T — oo with t/T finite, decaying with power law 
[19]. Berman 29|] has proved that in this regime the correlations are too weak to affect 
the MRHa distribution, which has the same FTG form as for a = 0. According to the 
simulation results shown in Fig. [H the MRHj distribution has the same Gaussian form (|32|) 
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throughout the interval < a < 1. Thus, the correlations appear to be equally irrelevant 
for the MRHi distribution. 

At a = 1 we enter a regime of stronger correlations. The correlation function {h{0)h{t)) 
diverges as In T at a = 1 and as T"~^ for a > 1. Although Berman's result on the irrelevance 
of correlations no longer holds at a = 1, the numerically determined MRHj distribution is 
still in good agreement with the Gaussian form (l32l) . 

Once we are well within the strongly correlated regime {a larger than ~ 1-2), the con- 
vergence with increasing N is much faster than for a < 1, just as in the case of the MRHa 
distribution pJi]. Our exact and numerical results MRHi for a > 1 are summarized in Fig. HI 
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FIG. 4: MRHi distribution for periodic boundary conditions and a in the range 1.2 < a < oo, 
scaled with the average, as in Eqs. (j28p and (j29p . For a = 2, 4 and oo the exact analytic results are 
shown, while the results for a = 1.2, 1.5, 3 and 8 were obtained from simulations with = 4096 
terms in the Fourier series ()24p . Note that the curves for a = 8 and a = oo are practically 
indistinguishable. 

A prominent feature in Fig. HI which we have already mentioned, is the increase in weight, 
with increasing a of the distribution function $i(a;) for small x, with a divergence at x = 
for a greater than a critical value ac- For small x, 

^i{x) ~ , X ^ , (62) 

with an exponent 7(a) that decreases monotonically as a increases, and changes from 
positive to negative at Oc- The Gaussian form of P{m,N) for < a < 1, with 



a= 1.2 
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(m)jv ~ cr]^^ ~ (InA^)-^/^, leads to $i(x) = 5{x — 1) on scaling by the average. This is 
consistent with 7(0;) = 00 for a < 1 in Eq. (1621) . The exact results in Eqs. ( 1371) . (H6l) . and 
(!60|) imply 7(2) = 1, 7(4) = — i, and 7(00) = — i. The curve for a = 3 in Fig. H] looks 
compatible with Oc = 3, but determining ac numerically with good precision requires a more 
extensive study, with careful attention to finite size effects near x = and a = 3. 

We have also studied the small x singularity for a = 5 and 6 numerically. In both cases 
the exponents 7 are in good agreement with the value — i, although the simulations can, 
of course, not exclude small deviations. These findings and the exact analytical results 
7(4) = — I and 7(00) = — i appear to indicate that 7(a) = — | for all a > au, where a„ is 
an upper critical value between 4 and 5. In the next Section we present some simple physical 
arguments which predict = 3 and = 5. We also obtain a formula for 7(a), a < 5 
which reproduces all the known exact values and agrees with the simulations. 

VI. CONNECTION BETWEEN THE MRH/ DISTRIBUTION AND THE DEN- 
SITY OF STATES NEAR EXTREMES 

As we saw at the end of Sec. II, for iid variables the density of states near extremes, 
studied by Sabhapandit and Majumdar 



20|, and the MRHj distribution are identical. Here 
we show that this equivalence continues to hold for correlated variables, as long as the 
probability distribution of the signals is time-translationally invariant. This is the case for 
correlated Gaussian signals with 1//° noise and periodic boundary conditions, and we use 
the equivalence with the density of states to study the small x singularity of 

In a time-translationally invariant system, all of the points traced out by the signal 
correspond to possible initial states. For a given realization of the signal, the distribution 
of heights measured from the maximum height is the same as the the distribution of the 
maximum height with respect to all these possible initial heights. Thus, averaging over 
all realizations yields identical distributions for the quantities /imax — h and /imax — h{0). 
We have checked the equivalence of the two distributions for periodic, correlated Gaussian 
signals both numerically and, for a = 2 and 4, by expressing both distributions in terms of 
propagators. 

We now present a simple picture that is very helpful in understanding the threshold be- 
havior of the density of states relative to the maximum. The idea is most readily understood 
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in the a — oo limit. In this case, considered in Sec. lIVDt each path is a smooth curve with 
a parabohc maximum, so that 5h = h{tm) — h(t) ~ (5^)^, where 6t = t — tm- Thus, the 
density of states relative to the maximum has the singularity p{5h) St/Sh for 
6h 0. Since the density of states is the same as the MRHi distribution, the exponent 7 
in Eq. fl62l) has the value — |, in agreement with the exact result in Eq. (1601) . 

As long as the path h{t) is smooth near its maximum, in the sense that typical maxima 
are parabolic, the above argument applies. Since the average increments of 1//" signals 
scale as 



for a > 1, the signals are expected to be twice differentiable for a > 5. Thus, we predict 
7 = — I for a > = 5, in agreement with the simulations for a =5 and 6. 

The lower critical value is the smallest value of a for which the 1//" signals are once 
differentiable, which, according to the scaling (l63l) is expected to be etc = 3. In this case, 
the paths near the maximum are basically rooftop- like, so that 6h ^ \6t\ and p ~ constant. 
This implies 7(8) = 0, which is compatible with the simulation results for a = 3 in Fig. |H 

To obtain a formula for 7(0;) in the interval 1 < a < 5, we assume that the scaling 
behavior in Eq. (163!) applies near the maximum of the paths. This implies p{Sh) 6t/6h ^ 
((5/i)(3-«)/(«-i) and so 7(a) = (3 - a)/{a - 1). 

In summary, the simple scaling picture predicts the small argument behavior 

p[bK) ~ bK< , $1(2;) ~ , (64) 

where 

a < 5 , 



7(a)=<|«-^ (65) 
"2 ' a > 5 , 

for both the density of near extreme states and the MRHj distribution. The above derivation 
may seem oversimplified, but the result is in remarkable agreement with those presented in 
Sees. IV and V. Expression fl65|) reproduces the exact values of 7(0;) for a =2, 4, and 00 
collected below Eq. (l62l) . Furthermore, the divergence 7 ^ 00 for a ^ 1 is in accordance 
with the fact that, in this limit, the MRHj distribution approaches a delta function centered 
at X = 1. It also agrees with all our simulations. We note, however, that a convincing 
numerical confirmation of 7(a) was not achieved in the regime 1 < a < 1.6 due to finite-size 
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corrections. We suspect the relation (l65i) is exact, but more effort is needed to put it on a 
firmer foundation. 

Tlie equivalence of the MRHj distribution and the density of near-extreme states and 
the scaling picture presented above furnish a simple physical explanation for the striking 
increase of both distributions, for small relative heights, with increasing a. As a increases, 
high frequency Fourier components of the Gaussian signal are suppressed. Typical signals 
become smoother and approach the maximum height less steeply, remaining close to the 
maximum for a longer time. Beyond a = 3 the approach to the maximum is basically 
tangential, becoming parabolic for a > 5. Thus, the density of near maximum heights 
6t/6h increases as a increases and diverges at the maximum for a > 3. A large density of 
near-maximum heights is equivalent to a high probability that the maximum height is close 
to the initial height. The scaling picture provides both a simple qualitative explanation of 
the singular small-argument behavior Eq. (!64|) of the two equivalent distributions p{Sh) and 
and a quantitatively accurate prediction for the exponent 7(a). 



VII. CONCLUDING REMARKS 



In this paper we have shown the importance of the reference point from which the max- 
imum is measured in extreme statistics. We found that the distributions of the maximum 
height relative to the average height and relative to the initial height are generally not the 
same. One reason for this is that both the average height and the initial height are fluctuat- 
ing quantities, but the distributions of their fluctuations are different. Furthermore, the two 
different reference heights impose different constraints on the paths consistent with a partic- 
ular value of the maximum height. For example, the probability for small maximum height 
relative to the average height is small, since only paths which remain close to the average 
for the duration of the signal contribute. In the case of a small height relative to the initial 
height, a much larger family of paths, which remain below the initial height most of the 
time, is allowed. This difference is highlighted by the analytical result that the MRHa dis- 
iribution has an essential singularity at small heights with exponentially suppressed values 



191], whereas the MRHj distribution has power law behavior near the origin. 



A notable feature of the aforementioned singularity of the MRHj distribution is that the 
exponent in the power law monotonically decreases with a, thus increasing the weight near 
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zero heights. This is in agreement with known persistence properties of 1//" processes [30| , 
where the weight of configurations persisting below the starting height shows a similar trend 
in a. To make this intuitively appealing connection more rigorous, further considerations 
would be required. 

The small height singularity of the MRHj distribution was described quantitatively with 
the help of a remarkable connection to the density of near extreme states. Namely, we found 
that they are identical for iid variables or, more generally, for signals drawn from time- 
translationally invariant distributions. It should be noted that we also considered boundary 
conditions other than periodic, with the aim of illustrating the dependence of boundary 
conditions on the extreme statistics. Since time-translational invariance is violated in these 
cases, the density of near extreme states remains open for further studies. 
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APPENDIX A: DERIVATION OF RESULTS FOR a = 4 

In the case a = 4, corresponding to random acceleration, the free space and half space 



propagators and Z\ are given by yM, 

Zoih,v;ho,Vo;T) = 3i/2(27r)-ir2 



X exp < —3t 



-3 



{h- ho- vt){h -ho- vot) + ^(^ - ^o)^^^ 



(Al) 



Zi{h,v;ho,Vo]T) = Zo{h,v] ho,Vo]T) + A{h,v; ho,Vo]T) , (A2) 

T,, , , I r dF r dG 

A{h,v;ho,Vo;s) = 



X 



271 Jo Fi/6 
exp [-{Fh + Gho) - (F-i + G'^ 



F + G 

X Ai {-F'/'v + F-^/^s) Ai {G'/^o + G~^/h) . (A3) 
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Here A(. . . ; s) denotes the Laplace transform £ A(. . . ; T) = dT e '^'^A(. . . ; T), and Ai{z) 
is the standard Airy function Ai 26|]. Using these results, we rewrite Eq. (|45|) in the form 

/oo 
dvA{Tn,v;m,v;s), (A4) 
-oo 

where indicates the inverse Laplace transform. 
It is convenient to calculate the moments 

d 



{771") = J dmm''^J^p,,{m,T) , (A5) 

which have simple Laplace transforms, and then construct the extreme distribution from 
the moments. After substituting Eqs. flA3l) and flA4l) into flASp . we first integrate over m 
and then over v, using 

/oo 
dvAi {-F'/h + F-^/h) Ai {G'/h + G-^/h) = 
-oo 



{F + G)"'/^Ai s{F-' + G-') 



-n2/3 



(A6) 



which follows from the integral representation (10.4.32) of Ai(2;) in Ref. 26j]. On expressing 
the right-hand side of Eq. ( ]A6p in terms of the Bessel function Ki/^iy) according to Eq. 



10.4.14) of Ref. 26|, making the changes of variables G = Fz, F = ^s^^'^y ^(1 



integrating over y and z with the help of Ref. [27j and Mathematica, and using C ^ [s °' ^] 
r{a + 1)~^T°', one obtains 



r(|z. + i) f T'/'Y _ r{lu + i) 

^ " ' ^ 18 r(z/ + l) 



r (1)-^ (m) 



r(z/ + i) 

which implies the scaled moments 

To calculate Pper(^, T), we consider the Laplace transform or generating function 



(A7) 



dme ^™Pper(m, 



n=Q 



(A8) 



Substituting Eq. (1A7I) into ( 1A8I) and carrying out some straightforward steps, one obtains 



r(i)"' we 
1 



n=0 

foo oo 

^ F n + 1 2 

n=0 ^ ' 

fOO 







(ixe-Vo (^2^r (I) ' (m)ex3/2 



(A9) 
(AlO) 

(All) 

(A12) 
(A13) 
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where Jo{x) is the standard Bessel function |26 |. 
Next we will need the relation 



[r'/' exp (-yr'/')] = m-'/'H{m^/'y) , (A14) 

where 

H{z) = y — \\ ^ , (A15) 

^r(n + l)r(|n + |) ' 

which follows from expanding the exponential function on the left-hand side of Eq. ( 1A14P in 

powers of y^'"^^^ and inverting the Laplace transform term by term. Evaluating the inverse 

Laplace transform of Eq. (]A12p using Eq. (]A14I) and changing the integration variable to 

z = m^^^y yields 

Pper{m,T) = m-^ dzH{z)Jo(^2^jT {ly^ m-^{m)z^/^^ . (A16) 



With the help of Mathematica one may express the series (lAlSP in terms of generalized 



hypergeometric functions j27| and evaluate the integral flA16p . This leads to 

H{z) = r(f)-\F4(|,ii,i;-Tfe-^) (Air) 

_ r(^\'^ r F f'2 2 7 4. 1_ 3\ 

^ \3) ^ 0-^4 1^3, 3, Q, 3, 108^ ) 

2 ^ 0-^4 I, -L) 35 2' 3' 108^ ) 

and our final result for the extreme distribution, 



Pper(m,T) = (m)-i<l>i(^) , (A18) 



where is given in Eq. (l46l) . 
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